Increasing of entanglement entropy from pure to random quantum critical chains. 
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It is known that the entropy of a block of spins of size L embedded in an infinite pure critical spin 
chain diverges as the logarithm of L with a prefactor fixed by the central charge of the corresponding 
conformal field theory. For a class of strongly random spin chains, it has been shown that the 
correspondent block entropy still remains universal and diverges logarithmically with an "effective" 
central charge. By computing the entanglement entropy for a family of models which includes the 
iV-states random Potts chain and the Zn clock model, we give some definitive answer to some recent 
conjectures about the behaviour of this effective central charge. In particular, we show that the ratio 
between the entanglement entropy in the pure and in the disordered system is model dependent and 
we provide a series of critical models where the entanglement entropy grows from the pure to the 
random case. 



The effects of quenched randomness on quantum crit- 
ical systems have been the subject of an intense experi- 
mental and theoretical activity for many years 1\. Dis- 
ordered one dimensional (ID) quantum models provide 
an important framework to study these effects because 
the quantum fluctuations are enhanced in this dimen- 
sion and the corresponding pure systems can be stud- 
ied using a wide variety of approaches, such as Bethe- 
ansatz solutions and conformal field theories techniques. 
For systems displaying strong spatial inhomogeneities, 
many results can be obtained by the application of the 
strong disorder renormalization group (SDRG) method 
introduced by Ma, Dasgupta and Hu For a 

wide variety of random quantum critical models, ranging 
from the random quantum chains with discrete symme- 
tries to the random Heisenberg models on fractal lattices 
H H, E S H ES [H E3, it has been shown that the 
effective disorder grows indefinitely at larger scales. This 
flow to strong disorder occurs in particular for the ran- 
dom spin 1/2-Heisenberg chain, the random XX chain 
and the random transverse field Ising model (RTFIM). 
The physics of these systems is captured by the so called 
infinite random fixed point (IRFP). Some of the main 
features of this point are a strong dynamical anisotropy 
and an extreme broad distribution of physical quantities 
which manifests through drastically different behaviour 
between average and typical correlation functions. 

Recently, Refael et Moore 13] have used the particular 
properties at the IRFP of these models to compute the 
entanglement of a block of spins of the corresponding 
chain. The entanglement of a segment of L spins with 
the remainder is defined as the von Neumann entropy of 
the reduced density matrix pl 



S(L) = —TrpL In p L . 



(1) 



In an infinite quantum critical chain without disorder, the 
entanglement of the segment with the remaining sites is 
proportional to InL with a coefficient depending on the 
central charge c of the associated conformal field theory 



EH03: 



S(L)~-lnL 



(2) 



The main result of |l3| is that the universal logarith- 
mic scaling of the entanglement entropy still holds at the 
IRFP of these models, where the conformal invariance is 
lost. In particular they showed that the entanglement 
entropy Srtfim, Srh and Srxx respectively for the 
RTFIM, the random Heisenberg and XX chain is: 

log 2 



Srh(L) = Sftxx(L) 
Srtifm(L) 
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Numerical evidences for the logarithmic scaling have 
been obtained for the random Heisenberg and XX chain 
in [HGjl. 

Analogously to the role played by the central charge, 
the prefactor of the logarithmic divergence can be a good 
estimator to characterise the universality class of the ran- 
dom chains. The continuum limit of the pure Heisenberg 
(H) and XX chains (XX) are described by the c = 1 
free boson theory (with different compactification radius) 
while the pure quantum Ising (TFIM) is described by 
the c = 1/2 free fermion theory From J5J), one has 
S H (L) = Sxx{L) ~ lnL/3 and S T fim{L) ~ lnL/6. 
Hence, the ratio between the random and pure logarith- 
mic prefactor is, for these two models, the same. This 
has suggested the intriguing possibility that the coeffi- 
cient appearing in the entropy scaling of any random 
fixed point is the product of the central charge of the 
associated pure critical system and an universal number 
(In 2). The findings ©, from which the entanglement 
of the pure systems turns out to be reduced by disor- 
der, suggest also a loss of entanglement along a general 
SDRG flow. Note that the loss of entanglement entropy 
was proposed in [2(J as a signal of the irreversibility of 
the renormalization group (RG) flow . This conjecture 
has been explored for the quantum Ising chain by vary- 
ing the transverse magnetic field and it was also found to 
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be consistent with the predictions for a Luttinger liquid 
with one single impurity |2l| . 

In this Letter we adopt the analysis of Refael and 
Moore to study the entanglement entropy of a family 
of random quantum critical chains, which includes the 
TV-state random quantum Potts chain and the Zjv clock 
model. By means of the SDRG procedure, we show that 
all these random chains are attracted by the infinite ran- 
dom fixed point of the RTFIM. Computing the associ- 
ated entanglement entropy, we show that that the cor- 
rection from a pure conformal fixed point to the random 
fixed point depends on N and thus is model-dependent. 
Moreover we find that for the random parafermionic Zn 
chains (defined later) with TV > 41, the entanglement is 
greater than the one in the pure system, thus providing 
a counter-example to the above conjecture. 

We consider a chain of spins with N states \si) = 
\0i), • • • , \(N — where i indicates the lattice site. 
The random quantum spin model we study is defined by 
the Hamiltonian H.n: 

jv-i 

i ri— 1 

N-l 

-J>£anr?. (4) 

i n—1 

in terms of the operators S z — e 2i ^' Nq S q ^qi (with q, q' = 
• • • N — 1), their hermitian conjugates S z and spin rais- 
ing operators T\s) = \s+l, mod A). The couplings Jj,i+i 
and the transverse fields hi are independent positive ran- 
dom variables drawn from some distribution. The co- 
efficients a n satisfy a n = a^-n to assure the Hamilto- 
nian to be hermitian and they are assumed to be disor- 
der independent. For each set of values {a n }, it exists 
a transformation from site variables to bond variables 



S S, 



i + 1 



S^'* which yields the same 



Hamiltonian with the J's and h's interchanged. The 
model @ is thus provided of a duality transformation. 
In the following we will always assume the couplings and 
fields to have initial equal distribution. 

If all the coefficients a n are equal to unity, a n = 1 
for n — 1 ■ ■ ■ N — I, the Hamiltonian (0J has a permuta- 
tional symmetry Sn and it defines the random quantum 
TV-states Potts model. The case a n = Sx,n corresponds 
to the random quantum Zn clock models. Senthil and 
Majumdar showed that these two models, with fer- 
romagnetic couplings, are attracted to the IRFP of the 
RTFIM. As they noticed, this implies that the scaling 
behaviour of many physical quantities (such as the mag- 
netisation or the mean correlation functions) depend on 
statistical properties entirely given by the IRFP distribu- 
tions which are N- independent. So, despite the fact that 
in the absence of disorder these systems have different 
type of correlations, the low-energy behaviour of the ran- 
dom systems presents ./V-independent scaling functions 



and exponents. On the other hand, the entanglement 
entropy at the IRFP of these models is, as we show later, 
sensitive to the number of spin states. This quantity can 
then discriminate between models with different N. 

Here we focus also on another class of random quan- 
tum spin chains defined by the Hamiltonian Q with 
a n = sin(7r/iV)/sm(7rn/iV). The interest in this model 
is motivated by the following facts: i) the corresponding 
pure model is critical at the self dual line J = h |23j and 
its fluctuations are governed by the Zn parafermionic 
field theories with central charge cjv = 2(N —l)/(N + 2) 
introduced in (the case N = 2 and N = 3 correspond 
to the Ising and the N = 3 state Potts model respec- 
tively), ii) The SDRG procedure can be performed for all 
this family of theories, iii) The scaling of entanglement 
entropy for the Zn spin chain can be computed using the 
properties of the IRFP and compared with its value at 
the parafermionic conformal critical points. Henceforth 
we refer to this model as the Zn parafermionic spin chain. 

We apply the SDRG approach to the model (@J by 
assuming the coefficients a n < 1. The renormalization 
equations are obtained by successively eliminating the 
maximum of the amplitudes of the bonds J^i+i and fields 
hi, getting an effective Hamiltonian for low energy de- 
grees of freedom. At each decimation step, if this maxi- 
mum turns out to be the field hi , the corresponding spin 

| Si) is frozen in the state |sj) = l/VFflOi) + |I t ) H h 

\(N — l)j)) which is the ground state of the dominant 
term — hi J2n=i Q «r". This generates an effective new 
coupling Jj_i ! , + i between adjacent sites. The coefficients 
a n enter into the renormalization as well. Degenerate 
second order perturbation theory yields: 



Ji—X.i 



Ji-lJi+1 ~ 2 K l (r\ 

+i = T ot n = a n — . (5) 

K\tli K n 



where the a n are the renormalizcd coefficients and 



i x- 1 



2nmn 
N 



(6) 



On the other hand, if the maximum is the coupling Ji,i+i, 
we set spin |sj) and spin |si+i) to have the same value, 
\si) = \si+i). The effective field h{ acting on the ferro- 
magnetic cluster |sj) ® |sj) is: 



hjh 



i'H+i 



KiJ, 



(7) 



as it could be guessed by duality. 

Due to our initial assumptions on the coefficients {a„}, 
we have k\ > 1 for N > 2. This assures that, during the 
renormalization procedure, the energy scale is decreas- 
ing: the fixed point of the problem is thus expected to be 
strongly attractive, even for weak disorder. It is impor- 
tant also to remark that the disorder drives the system 
away from the parafermionic critical point but the self- 
duality is preserved during the decimation procedure. 



J i-l,i+l 

FIG. 1: Schematic picture of the two kind of decimation step: 
a) Field decimation, b) Coupling decimation. 

Starting with the Potts values a n = 1 for n — ■ ■ ■ N— 1, 
the coefficients do not renormalize during the decimation 
procedure. In other words the permutational symmetry 
Sn is not broken along the flow of the couplings. On the 
contrary, by applying iteratively Eq0 and Eq. El the co- 
efficients {a n } defining the Zn parafermionic spin chain 
renormalize quickly to the values a n — ► 5\ tn . We can thus 
conclude that the random parafermionic Zn spin chain 
and the random Zn clock model share the same low- 
energy behaviour. In particular the couplings and fields 
distribution is attracted towards the IRFP distribution. 

The analysis of Refael and Moore can then be directly 
generalised to the more general random quantum Hamil- 
tonian The critical fixed point distribution yields 

a concentration nr of cluster at the log-energy scale T 
which scales as nr ~ l/T 2 . Hence the ground states of 
all these models present clusters which form over an av- 
erage length A ~ T 2 . Let us consider in detail a cluster of 
s spins. We want to compute the entanglement entropy 
of a subsystem of r spins with the remaining s — r spins 
of the cluster. If the cluster has been decimated by the 
field, it is frozen in the state \tl>): 

JV _ 1 s times 

M = -f=£l?>®l9>®"-®l?>- ( 8 ) 

* <7=0 

The entries of the reduced density matrix p r can be de- 
fined in the spin basis (j) % r , i,j = l..2 r as: 

(4\Pr\4) = E (41 ® {<A-Mm<t>\-r) ® Wr), (9) 
k,l—l 

where one traces over the basis (f>\_ r of the other s — r 
spins of the cluster. Using the wave function JSJl in 
one gets that the entanglement of a subsystem of r spins 
with the remaining s — r ones of a frozen cluster is log N 
for each r and s. Returning to the entanglement of a 
segment L embedded in an infinite chain, the frozen clus- 
ters which are totally inside or outside the segment L do 
not affect the total entanglement. The only contribution 
comes from the number of decimated clusters which cross 
the two edges, each of such clusters adding log TV to the 



FIG. 2: a) The state describing two spins decimated by the 
couplings, b) The state of a cluster of two spins after the 
field decimation, c) The contribution of the entanglement is 
given by the field decimated cluster which cross an edge of 
the segment (dashed line). In this example, we have three of 
such clusters. 

total value of the entanglement (see FigEJ. One has then 

S%° (L)~ 2 log Np{D) L (10) 

where S^ S (L) is the entanglement of a segment L in the 
disordered quantum chain (J3J , (D) l is the average num- 
ber of decimation which occur over an edge of the seg- 
ment L, p is the probability to have a field decimation 
instead of a coupling one. Due to the self-duality of these 
models, preserved along the SDRG, one has p = 1/2. 
The value of (D)l does not depend on N: it is entirely 
given from the fixed point distribution and it has been 
computed in 01 where it was found (D)l = 1/3 log L. 
Collecting all these results we finally arrive at the value 
of entanglement for the family of random models (0J with 
a„ < 1 : 

■ - log TV 

S%*(L)~-^-logL. (11) 

The entanglement <S^ S (L) depends thus on the num- 
ber of states N. Hence, differently from other quantities 
which are N- independent, it represents a good estima- 
tor to characterise the universality class of the random 
chains. Note that this quantity alone is not sufficient to 
identify a random phase transitions as it can be seen, for 
instance, from the fact that the random quantum Potts 
model with symmetry Sn have the same entanglement 
of the models with a lower symmetry Zn- This is anal- 
ogous to what happens in conformal field theory where 
different models can have the same central charge. 

Taking into account the value of the central charge cn, 
the entanglement S N ure (L) of the pure Zn parafermionic 
spin chain is 

nmre , . 2(N- 1) , x 

Comparing the values S N Llre and Sf} s (L), we can rule out 
the existence of a model-independent correction which 
links the values of the entanglement of the pure and of 
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the disordered system. Actually, there is no any clear 
relationship between the value of the entanglement pref- 
actor in the disordered chain and the central charge asso- 
ciated to the pure chain. It would be interesting to check 
other observables, for instance the scaling of the free en- 
ergy at the IRFP, which could clarify better if the value 
of this prefactor behaves as an "effective central charge" , 
as recently proposed 0, 0] . 
It is interesting to note that : 



foriV > 41, 



(13) 



thus invalidating the recent conjecture of a loss of entropy 
from pure to disordered critical chain. This conjecture 
followed along the lines of the c-theorem jiij which states 
the existence of a function of the couplings monotonically 
decreasing along the renormalization group. For this rea- 
son it is interesting to mention the fact that, in the pres- 
ence of disorder, the unitarity (one of the key hypotheses 
of the c-theorem) is in general lost. This is known to 
happen, for instance, in the case of the two-dimensional 
random bond 3-states Potts model which flows towards 
a random fix ed p oint where the conformal symmetry is 
restored [2ol l26| . In this case, one can show that the 
energy-energy operator, which determines the properties 
of the RG equations, gets a negative norm in the disor- 
der limit. This explains the fact that, for this model, the 
central charge associated to the disordered random fixed 
point is greater than the one associated to the pure fixed 
point. On the other hand, the loss of unitarity does not 
necessarily imply the increasing of the central charge, as 
it has been explicitly shown for instance in |27j . where a 
particular class of disordered conformal field theories was 
studied. In an analogous way, we think that, in the pres- 
ence of disorder, one cannot predict the behaviour of the 
entanglement on the basis of a generalised c— theorem, as 
our results clearly show. 

To conclude, in this Letter we have shown that the en- 
tanglement entropy of a class of random quantum criti- 
cal chains, including the random quantum Potts and the 
random quantum clock chain, obey an universal loga- 
rithmic scaling. We have confirmed thus the entangle- 
ment as an important quantity to characterise (although 
not completely) the universality class of this entire se- 
ries of disordered critical models. Moreover our results 
provide a definitive answer to some recent hypotheses on 
the entanglement entropy in the random quantum critical 
chains. By considering the Zn parafermionic spin chains, 
we could compare the value of the central charge of the 
pure system with the value of the entanglement entropy 
of the disordered one for each value of N. We found that 
the ratio between the central charge and the prefactor 
of the entanglement logarithmic scaling of the disordered 
models depends on N. We could thus exclude the exis- 
tence of a model-independent correction between these 
two quantities. Moreover, we provide some explicit mod- 
els where the entanglement of the pure system is smaller 



than the entanglement of the disordered one. Thus the 
recent conjecture about a generalised c-theorem concern- 
ing the entanglement does not hold when the disorder is 
present. 

I would like to thank Daniel Cabra, Mirko Cinchetti, 
Vladimir Dotsenko, Andreas Honecker, Nicolas Lafloren- 
cie, Pierre Pujol and Alberto Rosso for stimulating and 
interesting discussions. I also thank the CNRS for finan- 
cial support. 



[1 

[2 

[3. 

[4 

[5 
[6 

[-: 

[8 

[9 

[10; 

[11 

[12 

[13 

[14 

[IS 

[16 

[17 
[18 
[19 

[20 

[21 
[22 

[23 
[24 
[25 

[26 
[27 



Majumdar, Pliys. Rev. Lett. 76, 



Electronic address: santa@lptl.u-strasbg.fr 
For a review, see: R. N. Bhatt, in Spin Glasses and 
Random Fields, edited by A. P. Young (World Scientific, 
Singapore, 1998) p. 225. 

For a review, sec: F. Igloi and C. Monthus, 
cond-mat/0502448 

S. K. Ma, C. Dasgupta and C. K. Hu, Phys. Rev. Lett. 
43, 1434 (1979); C. Dasgupta and S. K. Ma, Phys. Rev. 
B 22, 1305 (1980). 

D. S. Fisher, Phys. Rev. B 22, 1305 (1980), Phys. Rev. 
Lett. 69, 534 (1992); 

D. S. Fisher, Phys. Rev. B. 50, 3799 (1994) 

R. A. Hyman and K. Yang, Phys. Rev. Lett. 78, 1783 

(1997). 

C. Monthus, O. Golinelli and Th. Jolicoer, Phys. Rev. 
Lett. 79, 3254 (1997). 

G. Refael, S. Kehrein and D. S. Fisher, Phys. Rev. B. 66, 
060402 (2002). 

J. .1. Latorre, R. Orus, E. Rico, and J. Vidal, Phys. Rev. 
A 71 064101 (2005). 

D. Carpentier, P. Pujol, and K. U. Giering, 
cond-mat/0506597 
T. Senthil and S. N. 
3001 (1996) 

R. Melin, B. Doucot and F. Igloi, Phys. Rev. B. 72, 
024205 (2005) 

G. Refael and J. E. Moore, Phys. Rev. Lett. 93, 260602 
(2004). 

A. Galindo, M. A. Martin-Delgado Rev.Mod.Phys. 74 
(2002) 347-423; quant-ph/0112105 

C. Holzjey, F. Larsen and F. Wilczek, Nucl. Phys. B 424, 
443 (1994) 

G. Vidal, J. I. Latorre, E. Rico and A. Kitaev, Phys. 
Rev. Lett. 90, 227902 (2003) 

P. Calabrese and J. Cardy, J. Stat. Mech. 06, 002 (2004). 
N. Laflorencie, Phys. Rev. B. 72, R140408 (2005). 
G. De Chiara, S. Montenegro, P. Calabrese and R. Fazio, 
cond-mat/0512586 . 

J. .1. Latorre, C. A. Lutken, E. Rico, and G. Vidal, 

Phys.Rev. A 71 034301 (2005). 

G. Levine, Phys. Rev. Lett. 94, 062301 (2005). 

V. A. Fateev and A. B. Zamolodchikov, Sov. Phys. JETP 
62, 215 (1985) 

F. C. Alcaraz, J.Phys.A:Math.Gen. 20, L623 (1987). 
A. Zamolodchikov, JETP Lett. 43, 730 (1986). 

VI. Dotsenko, M. Picco and P. Pujol, Nucl. Phys. B455, 
701 (1995). 

J. L. Jacobsen and M. Picco, Phys. Rev E 61 R13 (2000) 
VI. Dotsenko, X.S. Nguyen, and R. Santachiara, Nucl. 



Phys. B 613, 445-471 (2001). 



